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Abstract 

A scheme for generating Schrodinger cat-hke states of a single-mode optical 
field by means of conditional measurement is proposed. Feeding into a beam 
splitter a squeezed vacuum and counting the photons in one of the output 
channels, the conditional states in the other output channel exhibit a number 
of properties that are very similar to those of superpositions of two coherent 
states with opposite phases. We present analytical and numerical results 
for the photon-number and quadrature-component distributions of the con- 
ditional states and their Wigner and Husimi functions. Further, we discuss 
the effect of realistic photocounting on the states. 



1 Introduction 



It is well known that according to the basic-theoretical principles of quantum- 
mechanics the superposition of macroscopically distinguishable quantum states 
can give rise to quantum interferences such that the resulting states are highly 
nonclassical. Schrodinger [1] illustrated this phenomenon by a gedankenex- 
periment to get a cat into a superposition of a live and a dead cat. A number 
of systems have been studied with the aim of the realization of Schrodinger 
cat-like states, where the "cat" is typically a mesoscopic system that has both 
microscopic and macroscopic (i.e., classically distinguishable) properties. 

In harmonic oscillators typical examples of Schrodinger cat-like states 
are superpositions of two coherent (i.e., most classical) states with opposite 
phases [2] . The superposition states exhibit some properties similar to those 
of simple statistical mixtures, but they also reveal typical interference fea- 
tures. Measuring the quadrature-component distribution, one observes two 
peaks that change their mutual distance in dependence on the phase of the 
quadrature component until they eventually overlap. In this particular case 
the difference between a coherent superposition and a statistical mixture is 
the most distinct. Whereas in the former case quantum interferences are ob- 
served which give rise to an oscillatory behavior of the quadrature-component 
distribution, in the latter case a single peak without interference structure is 
observed. The creation of such states in realistic experiments is not trivial. 
Proposals for preparing vibrations in molecules or crystals in Schrodinger 
cat-like states have been made [3, 4]. Recently, it has been proposed that 
Schrodinger cat-like states of a single harmonically bound (trapped) atom 
can be produced by appropriately driving the atom [5], and experiments 
have successfully been performed [6] . The possibility of preparing a harmon- 
ically bound atom in a superposition of two coherent squeezed states with 
opposite phases has also been studied [7]. 

Several proposals have been made to prepare a single-mode radiation field 
(which also corresponds to a harmonic oscillator) in a Schrodinger cat-like 
state (see, e.g., [8] and references therein). In particular, it has been proposed 
that conditional measurements may advantageously be used to realize such 
states. It is well known that if a quantity of a subsystem of a correlated 
two-part system prepared in some entangled state is measured, the state of 
the other subsystem "collapses" to a particular state. To produce conditional 
states of the type of Schrodinger cat-hke states, the use of a scheme for optical 
back-action-evading measurement in nonlinear media [9] was suggested [10]. 
The calculations show that when the photon number of the readout mode is 
measured, then a superposition of macroscopically distinguishable quantum 
states is generated in the signal mode [10]. To improve the scheme, it was 
proposed that a squeezed vacuum at the signal frequency is injected instead of 
amplifying the signal after back-action-evading measurement [11]. Recently, 
a modification of this scheme was studied, with special emphasis on the 
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experimental feasibilities with current technologies [12, 13]. 

In this paper we show that Schrodinger cat-like states can already be 
obtained using a simple beam-splitter scheme for a conditional measurement 
of the type considered recently in [14, 15]. The calculations show that when 
a squeezed vacuum is injected in one of the input channels (the second input 
channel being unused) and the photon number of the mode in one of the 
output channels is measured, then the mode in the other output channel is 
prepared in a conditional state that has the typical features of a Schrodinger 
cat-like state. In particular, the conditional states can be regarded as super- 
positions of two quantum states that are well localized in the phase space 
and bear a strong resemblance to squeezed coherent states. To demonstrate 
this, we analyse the states in terms of the photon-number and quadrature- 
component distributions and the Wigner and Husimi functions. We further 
discuss possible modifications of the states in realistic multichannel photon 
detection. 

The paper is organized as follows. In Sec. 2 the conditional states are 
derived, and in Sec. 3 their properties are discussed. Section 4 is devoted 
to a realistic detection scheme. Finally, a summary and some concluding 
remarks are given in Sec. 5. 



2 Basis equations 

Let us consider a lossless beam splitter and assume that the input fields can be 
regarded as being effectively single mode fields, with photon destruction and 
creation operators and a],, respectively {k = 1, 2). The photon destruction 
and creation operators of the output modes bk and bl, respectively, can then 
be obtained using the well-known input-output relations 
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bk — ^ Tk,k' O'k'i (1) 
fe'=i 

where 

.j^J cosOe^^- sin^e^^« A 

is a SU(2) matrix whose elements are given by the complex transmittances 
and refiectances of the beam splitter from the two sides. Equation (1) corre- 
sponds to a unitary transformation of the operators in the Heisenberg picture, 
hk = VakV^ {k = 1,2). Equivalently, the Schrodinger picture can be used, in 
which the photonic operators are left unchanged and the density operator 
is transformed such that the output-state density operator ^out is obtained 
from the input-state density operator as 



^'out 



= V^'QinV. (3) 
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The operator V can be given by [16, 17] 
where 

L2 = l{a\a2 — alai) , Ls = \{a\ai — 0,20,2) , (5) 

and v^o = (note that ipo is a global phase factor that may be omitted without 
loss of generality) . 

Now, let us assume that the mode in the first input channel is prepared in 
a state described by a density operator ^ini and the second channel is unused, 
so that the input-state density operator reads as 

^in = ftni ® |vac2)(vac2|. (6) 

Using Eqs. (4) and (6), the output-state density operator (2) can be given 
by [14, 15] 



^out = 13 13 { 

712=0 m2=0 



y/m,2\n2\ 



m2+n2 



r|2'r^"'ftni|T|"^"Hal)"^e-^'^^"i"i ® |m2)(n2| 



(7) 

where |T| = cos^ and |i?| = sin^, and \nk) are the eigenstatcs of the photon- 
number operators a^Cfc. From Eq. (7) we see that the output modes are, 
in general, highly correlated. When the photon number of the mode in the 
second output channel is measured and m2 photons are detected, then the 
mode in the first output channel is prepared in a quantum state whose density 
operator ^outi reads as 

^outl ^2) = ;^r77 p , TT- 8) 

Tri((m2|^out|"22)) 
The probability of such an event is given by 

P{m2) = Tri((m2|^out|"^2)) 

= E (""'Vl-lTlT'lTr^— )(ni|^iniK). (9) 

In particular, if the first input mode is prepared in a squeezed vacuum 
state, we may write 

ftni = ^(0|vaci)(vaci|5t(0, (10) 
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where 



^(Olvaci) = exp{-i [{^a\f - Ta?]} |vaci 



ni=0 



2"i ni! 



(11) 



^=|^|e'"^«, K = e"^« tanh |^|. Combining Eqs. (7) and (8) and using Eqs. (6), 
(10), and (11), we derive that 



where 



^outl("^2) = |*r„2)(*m2|, 
I ^"12) |^m2('^)) 



-j^ 00 

-7=^ E Cm2,ni(a) 
\j-Nm2 "1=0 



(12) 
(13) 



(a) 



(ni + 7722)! 



i(ni + 7712) + 1 



X I [1 + (-l)"i+™2] (la) 



'77i! 



(14) 



a = laje*'''", with |a| = |T|2|k| and ip^ = '^^t + V'g- In what follows we will 
restrict attention to real values of a (— 1 < « < 1), i.e., = 0, tt, since from 
Eqs. (13) and (14) the effect of other phases ip^ is simply a rotation in phase 
space. Applying the relations (A 2), (A 3), and (A 6) in the appendix, the 
normalization constant 

00 

|2 



■^m2 ^ ] |Cm2,' 
ni=0 



,ni I 



(15) 



can be given by 



m2 



(1- 

■.m2\ 



m2 



(7712!) 



£ (m2-2A;)!(A:!)^(2«)2^' 



(16) 



where the symbol [x] in the summation upper limit denotes the integral part 
of X. Similarly, using Eqs. (6), (10), and (11) and applying the relations (A 
2), (A 3), and (A 6), Eq. (9) yields 



P(7772) 







(l-|Tp) 




\T 


2(1 -q;2) 



m2 



jm2 



7772! 



£ (m2-2A;)!(A;!)2(2Q;)2'^- 



(17) 



5 



From Eqs. (13) and (14) we easily see that when the detected number 
of photons in one output channel, 1712, is even (odd), then the mode in the 
other output channel is prepared in a quantum state |^m2) that contains 
only contributions from photon-number states with even (odd) numbers of 
photons. This property, which gives rise to oscillations in the photon-number 
distribution of the output state \'$„i2): obviously reflects the fact that the 
squeezed vacuum that is fed in consists of pairs of photons. In particular, 
when the number of detected photons, m2, is zero, then the output state is 
again a squeezed vacuum, but with the parameter a in place of k. 

3 Properties of the conditional states 

To study the properties of the conditional states = |\E'm(ci)) in more 

detail, we will calculate the photon-number and quadrature-component dis- 
tributions and the Wigner and Husimi functions. Further, we will show that 
the states |^m) can be represented as superpositions of two macroscopically 
distinguishable "quasicohcrcnt" squeezed states. For notational convenience 
we will omit the subscripts 1 and 2 introduced above to distinguish between 
the two output channels. 

3.1 Phot on- number distribution 

Recalling Eq. (13), the photon-number distribution 



where Cm,n{(^) and Mm are given in Eqs. (14) and (15), respectively. In 
particular, the mean photon number 




(18) 



of a state l^'m) reads as 




(19) 



00 




(20) 



n=0 



can be given by 



2 ^ kak^m X] ^'^r 



-I -1 



1 — a 



- + m 




k=0 k=0 



(21) 
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where ak^m = {2a)^'^^/[{m — 2k)\ (kl)"^]. Examples arc shown in Fig. 1(a). 
We see that the number of photons that can be found in increases with 
m. This is simply a consequence of the beam splitter transformation. Since 
one of the input channels is unused, the mean numbers of photons in the 
two output channels are proportional to each other, the ratio being given by 
|T//?p. Note that when no photons arc detected. m = 0, then (n) reduces to 
the mean number of photons of a squeezed vacuum, (n) = q;^/(1 — a^). 

A measure of the deviation of the photon-number distribution from a 
Poissonian is the Mandel Q parameter [18] 



(n) 



1, 



which can be given by 



a" d 



^ {n) da^ 



log 



a' 



1. 



(22) 



(23) 



The dependence on a and m of Q is shown in Fig. 1(b). In particular, for 
even m we find that Q > for all values of a, which means that the photon- 
number statistics of \^m) is super- Poissonian. For odd m and small values 
of |a| the statistics become sub-Poissonian ((5<0). Note that the behavior 
is typical for Schrodinger cat-like states [8]. 



3.2 Quadrature distributions 



In order to calculate the conditional quadrature-component distribution (i.e., 
the phase-parametrized field-strength distribution) 

p(x,(^|m) = (24) 

which can be measured in balanced homodyne detection, we first expand the 
eigenvectors jx, </?) of the quadrature component 



x{ip) = 2-5 {e-'^^a + e'^^a^) 
in the photon-number basis as [19] 



\x. 



^) = (.)-.exp(-i.^)E^H„ 



{x)\n) 



(25) 



(26) 



(H„ is the Hermite polynomial). Using Eqs. (13) and (26) and applying the 
relations (A 3) and (A 6), the conditional quadrature-component distribution 
(24) reads as 



p(x, ip\m) 



\a\ 



■exp 



■ X 



X 



(27) 
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where the abbreviation 

A = 1 + - 2q; cos(2(/?) (28) 

has been used. 

Prom Fig. 2 we see that for (/? near 7r/2 the quadrature-component distri- 
bution p{x, Lp\m) (m> 0) exhibits two separated peaks, whereas for close to 
or TT an interference pattern is observed. It should be noted that although 
the quadrature-component distribution (27) bears a strong resemblance to 
that obtained in Ref. [11], a more detailed comparison shows that they are 
different. Clearly, the beam splitter transformation, Eqs. (3) - (5), can not 
be identified, in general, with the transformation in the back-action-evading 
scheme considered in Ref. [11]. 



3.3 Wigner function 

Using Eqs. (13) and (26) together with the relations (A 3) and (A 6), the 
Wigner function W{x,p\m) of the state l^'m). 



W{x,p\m) = - [ dye^'Py{x-y\'i/m){^m\x+y), 

TT J 

— oo 

can be calculated in a straightforward way. We obtain 



W{x,p\m) 



7r3/W„.[2(o; + l)] 



2 +00 ^ 

— oo ^ 



a 



1+a 



[x+y) 



1 — a 
1 + a' 



which after calculation of the integral [20] can be given by 



W{x,p\m] 



exp ( —Xx^ — ^ 



X 



E 



where 



[im] 



H 



■m—k 



i\/aX i X + i 



(2|«|) 



m—2k 



(29) 



(30) 
(31) 



(32) 



(33) 



to im-2k)\{k\f 

Equation (32) reveals that when no photons are detected, m = 0, then the 
Wigner function exhibits a single peak and simply corresponds to a squeezed- 
vacuum Gaussian, as already mentioned in Sec. 2. Examples of W{x,p\m) 
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(m > 0) are shown in Fig. 3. From the plots two separated peaks and an 
oscillatory regime between them can be seen - a shape that is typical of 
Schrodinger cat-hke states. The separation of the two peaks is seen to in- 
crease with the number m of detected photons. Since with increasing m the 
number of photons in the conditional state |\E'm) also increases, which is a 
consequence of the beam splitter transformation, the behavior is quite similar 
to that of a superposition of two coherent states. 



3.4 Husimi function 

The Husimi function Q{x,p\m) of the state l^'m) is defined by 

Q{x,p\m) = j-m^m)\', (34) 

where is a coherent state, and /3 = 2~^/^(x + ip). It is worth noting 
that, contrary to the Wigner function, the Husimi function is a phase-space 
function that can directly be measured in multiport balanced homodyning 
using six-port [21, 22] or eight-port schemes [23, 24]. Expanding \(3) in the 
Fock basis. 



oo 



'^'E^ln), (35) 



=0 Vn^- 

and recaUing Eq. (13) and the relations (A 3) and (A 7), the scalar product 
can easily be calculated. After some algebra we find that Q{x,p\m) 
can be written as 

Q{x,p\m) = ^^^^ \Y{„,\j^i^{x + ip)\ 

X exp{-i[(l-a;)x2 + (l + Q;)p2]}. (35) 

As expected, for m = the Husimi function is Gaussian, whereas for m > 
a two peak-structure is observed. Since the Husimi function is always non- 
negative, the interference properties of the state are not so apparent as in 
the case of the Wigner function. 



3.5 Component states 

Schrodinger cat-like states are commonly defined as superpositions of two 
macroscopically distinguishable quantum states. From Eqs. (13) and (14) it 
is seen that \^m) can be regarded as a superposition of states l^l'^'') and 
as follows: 

+ (37) 
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where 



with 



M^^ n=0 

(n + m)! 



WW r[(n + m)/2 + l]V^^^2« 



n+m 



(38) 



(39) 



The normahzation factor is calculated to be [see Eqs. (B 3) and (B 11)] 



v^2™r(m + 3/2) 



xF 



i(m + l), |(m+l), i(2m + 3), 1 - a' 



(40) 



A/"™ being given in Eq. (16). In Eq. (40) F(a,/5,7, 2;) is the hypergeometric 
function, and the normalization constant A in Eq. (37) is simply given by 
A = \{M^yUmY''^. Note that ^2Mm, i.e., for larger m (the 

approximation is very good for m>4). 

Plots of the Wigner function of the states are given in Fig. 4. 

The behavior of the states |^^') is quite similar. From the figure we see 
that the states l^^-*) (and also the states 1^^'')) ^^"^ very close to squeezed 
coherent states. For chosen a the squeezing effect decreases with increasing 
m. The small deviation from Gaussian states is indicated by the negative 
values of the Wigner function (which in our case are of the order of magnitude 
of — 10~^). To illustrate the difference between the states l^l/^-*) and the 
squeezed coherent states, let us consider the Husimi function Q^'^\x^p\m) of 
1^^)). According to Eq. (B 13) we obtain 



>{x,p\m) = , . e I''! 



X exp 



ia(/32 + /3*')l D_^_i(±V^/3) 



(41) 



where (3 = {x + ip) / and 0^(2;) is the parabolic cylinder function (which 
indicates the deviation from a Gaussian). Recalling the asymptotic behavior 
of the D_m_i{z) for large values of m [26], we find that 



Q^'^\x,p\m) 

oc exp 



■|/?|^ ± y/omi((3 + (3*) - \y/omi\ 



(42) 



(m — >• 00), which corresponds to the scalar product |(/3|7)P between the 
coherent states and I7), with 7 = ±y/am. Thus for sufficiently large 
numbers of detected photons, m, the component states |^^^) approach co- 
herent states. This kind of behavior has also been found in Ref. [11] for the 
states studied therein. 
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4 Realistic photon counting 



In order to produce the conditional states l^m), highly efficient and precise 

photocounting is needed. Unfortunately, there have been no highly efficient 
photodetectors available which precisely distinguish between m and m+1 pho- 
tons. Recently the proposal has been made to measure the photon-number 
statistics using photon chopping [25]. The mode to be detected is used as 
an input of an optical 2A^-port, the other input ports being unused. To 
each of the output ports a highly efficient avalanche photodiodc is connected 
which can distinguish between photons being present or absent. The photon- 
number statistics of the input mode can then be obtained from the recorded 
output coincident-event statistics. 

In particular, if m photons are present in the input, the probability of 
recording k coincident events is given by [25] 



(43) 



for k < m, and P/v(A;|m) = for k > m. Note that P/v(A;|m) — > 6k,m for 
N^oo. In Eq. (43) perfect detection is assumed. The effect of nonperfect 
detection corresponds to a random process such that photons are excluded 
from detection with probability 1^, rj being the efficiency of the photodiodes 
(0 < 77 < 1). The probability of recording k coincident events then modifies 
to 

PNAk\m) ^J2Mk\l) Mi,m{v): (44) 
I 

where the matrix Mi^m^rj) is given by 



MUv)=[jy{l-Vr~^ (45) 

for l<m, and Mi^m{r})—0 for l>m. 

Since detection of k coincident events can result from various numbers m 
of photons, the conditional state is in general a statistical mixture. Therefore 
in place of Eq. (12) we now have 

Qoutik) ^Y.PN,r,(m\k) (46) 

m 

where l^m) is given in Eq. (13), and Pjv.T,(m|A;) is the probability of m photons 
being present under the condition that k coincidences are recorded. The 
conditional probability PN,r]{'>Ti\k) can be obtained using the Bayes rule as 

PN,M\k) = -^^pNMm) Pirn). (47) 
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Here P{m) is the prior probability (17) of m photons being present, and 
accordingly, P^^rii^) is the prior probability of recording k coincident events, 

PN,r,{k) ^^pN,r,{k\m) P{m). (48) 

m 

Examples of the Wigner function and quadrature-component distribu- 
tions of the mixed state (46) are plotted in Figs. 5 and 6, respectively. We 
see that the quantum interferences can still be preserved also for reahstic val- 
ues of the number of photodiodes and their efficiencies, such as A'" = 10 and 
?7 = 80%. The interference structure is (for chosen and rf) more and more 
smeared with increasing k. Clearly, a larger number of detected coincidences 
implies (for chosen and r]) a larger probability of "losing" some of the 
photons. Any lost photon switches the parity of the conditional state (from 
even to odd and vice versa) and therefore destroys the interference pattern. 

5 Conclusion 

We have shown that Schrodinger cat-like states can be generated by condi- 
tional measurements using a simple beam splitter scheme. When a squeezed 
vacuum and an ordinary vacuum are mixed by a beam splitter and in one 
of the output channels the number of photons is measured, then the condi- 
tional quantum state in the other output channel reveals all the properties 
of a Schrodinger cat-like state. 

To demonstrate this, wc have analysed the conditional states in terms of 
the photon-number and quadrature-component distributions and the Wigner 
and Husimi functions and have presented both analytical and numerical re- 
sults. We have studied the component states and shown that they are very 
close to squeezed coherent states and approach coherent states for sufficiently 
large numbers of detected photons. We have found that the basic features of 
Schrodinger cat-like states, such as the appearance of two separated peaks 
and the interference pattern, become more pronounced for larger values of 
the transmittance of the beam splitter, the squeezing parameter of the input 
state, and the number of detected photons. On the other hand, increasing 
transmittance implies a decrease of the probability of photons being present, 
so that the "better" Schrodinger cat-like states appear more rarely. 

We have also discussed the problem of producing the Schrodinger cat-hke 
states under the conditions of realistic photocounting. For this purpose we 
have assumed multichannel detection using highly efficient avalanche pho- 
todiodes. As expected, the measurement smears the interference structure. 
However, for properly chosen parameters the interference structure can still 
be found even for a realistic arrangement of detectors. 
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Appendix A Sum. Rules 

Using the expansion in the Fock basis of |\&m)) Eq. (13), various photon- 
number summations must be performed in the further calculations, which 
can advantageously be done by means of Hermite polynomials. For this 
purpose we rewrite the coefficients Cm.2,niioi), Eq. (14), as 



Cm2,ni (<^) 



[2{n + k + 5)]\ 



{n + k + S)\J{2n + 6)\ 



n+k+S 



(A 1) 



where 6 = for ni = 2n and m2 = 2k. and 6=1 if ni = 2n+l and m2 = 2k+l. 
Note that Cm2,ni(Q') — otherwise. Recalling the relation 



H2„(0) = (-l)-(2n)!/n!, 



we see that 



-'m2,ni 



H 



2(n+k+S) 



(0) 



n+k+S 



(A 2) 



(A 3) 



This enables us to apply standard summation rules, such as Mehler's formula 



CO 1 ^ 



k=0 



1 



: exp 



2xyz — (,T^ + y'^)z'^ 



(A 4) 



Taking the Ith and jth derivatives with respect to x and y, respectively, of 
both sides of this equation and using the relation 



(A 6) 



we derive that 



^1 1 

E^H,+K^)H.+,(?/)(i^) - (i_^2)(/+,-+l)/2 



X exp 



(i-z^Y 



1-^2 

X — zy ^ 



H 



j-k 



y — zx 



(A 6) 
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Another useful sum rule is 



k=0 



(A 7) 



which may be derived by taking the jth derivative of the generating function 
of the Hermite polynomials, 



k=0 



(A 8) 



Appendix B Derivation of the relations (40) 

and (41) 



Prom Eqs. (13), (14) and (37) - (39) we see that 



Cm,n('-'^) — 2 



.(-) 



and hence 



E lcm,n(«)r = i E + E I41(«) 



n=0 



+ 2Re 



n=0 



E4U(«)4-i(«) 



.n=0 



Thus, 



•Ae^ = El41U«)^=2A4^-x„, 



n=0 



where is given in Eq. (16), and 



fB 11 



(B 2) 



(B 3) 



^m = E(-i 

n=0 



Combining Eqs. (B 4) and (39) yields 

THm + n + 1] 



for2[(m + n + 2)/2]n! 



Using the relation 



r(2n) 



r(n + l/2) 20F 



r(n). 



n+m 



(B 4) 



(B 5) 



(B 6) 
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from Eq. (B 5) we obtain 



n=0 



-2|q;| 



\n+m 



(B 7) 



Substituting in Eq. (B 7) for the F function the integral representation, we 
may write 



-2\a\y 



oo oo 



TT 



dti / dt2 







m-l)/2 



n=0 



n! 



which reduces to 



T = 



[-2\a 



TT >/ L 


oo 



We first calculate the t2 integral [26] to obtain 
m\{-^/2\a\)'^ 



T — 



TT 



Calculating the resulting ti integral [26] , we eventually obtain 

mP(-|a|)'" 



^"^ 0r2"»r(m + 3/2) 



X F 



|(m+l), i(m + l),m + |, 1 - 



(B 8) 



(B 9) 



(B 10) 



(Bll) 



where F{a, f3,'y, z) is the hypergeometric function. 

The calculation of (^^^1/9) can be performed in a similar way. Using 
Eqs. (13), (14), (35), and (B 6), we have 



(*L^^I/3) = (±1) 



oo 
n=0 



±V2a/3)" 
n! 



|(n+m+l) 



(B 12) 
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Substituting in Eq. (B 12) for the F function again the integral representation, 
we arrive at an integral of the type of the t2 integral in Eq. (B 9), which may 
be calculated using standard formulas [26]. We then obtain 




X e 



— m— 1 



(B 13) 



where D^{z) is the parabolic cylinder function. 
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(a) 



(h) 
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Figure 1: (a) Mean photon number (n) of the conditional state as a 

function of a for various numbers m of measured photons, (b) Dependence 
on a of the Mandel Q parameter for the same values of m as in (a). 
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Figure 2: Quadrature distribution p{x, (f\m) of the conditional state |^m) for 
a = 0.6 and various numbers m of measured photons [(a) m= 1, (b) m = 2, 
(c) m = 3, (d) m = 4]. 
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(c) (d) 

W[x,p\m) W{xMm) 




Figure 3: Wigner function W{x,p\m) of the conditional state \^m) for a = 0.6 
and various numbers m of measured photons [(a) m=l, (b) m = 2, (c) m = 3, 
(d) m = 4]. 
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Figure 4: Wigner function W^'^\x,p\m) of the state l^^-*) for a = 0.6 and 
various numbers m of measured photons [(a) m=l, (b) m = 3]. 
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(c) (d) 

Wix,p\k) W{x,p\k) 




Figure 5: Wigner function of the conditional state with reahstic photodetec- 

tion {N = 10, ri = 0.8) for k, = 0.75, |Tp = 0.8 {a = 0.6) and various numbers k 
of coincident events [(a) k = l, = 10.99%; (b) k = 2, pNjiik) = 2.95%; 

(c) k = 3, pN,-n{k) = 0.69%; (d) = 4, Piv,^(A;) = 0.16%], .the probabihties 
PN,ri{k) of the coincidences being calculated according to Eq. (48) 
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Figure 6: Quadrature distribution of the conditional state with reahstic pho- 
todetection (A^=10, 77 = 0.8) for the phase parameters (/? = (fuh hne) and 
(p = 7r/2 (broken hne), various numbers k of coincident events [(a) k = 2, (b) 
A; = 3], and k = 0.75, |Tp = 0.8 {a = 0.6). 
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